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Abstract 

We show that the initial value problem of a periodic box-ball system 
can be solved in an elementary way using simple combinatorial methods. 

A periodic box-ball system (PBBS) is a dynamical system of balls in an 
array of boxes with a periodic boundary condition 0[2j- The PBBS is obtained 
from the discrete KdV equation and the discrete Toda equation, both of which 
are known as typical integrable nonlinear discrete equations, through a limiting 
procedure called ultradiscretization [3|E]- Since the ultradiscretization preserves 
the main properties of the original discrete equations, and the solvability of the 
initial value problem being an important property of integrable equations, we 
expect that the initial value problem of the PBBS can also be solved. In fact, the 
initial value problem for the PBBS was first solved by inverse ultradiscretization 
combined with the method of inverse scattering transform of the discrete Toda 
equation and recently by the Bethe ansatz for an integrable lattice model with 
quantum group symmetry at the deformation parameter q = and g = 1 0. 
These two methods, however, require fairly specialized mathematical knowledge 
on algebraic curves or representation theory of quantum algebras. 

An important property which characterizes a state of the PBBS is the fun- 
damental cycle of the state, i.e., the length of the trajectory to which it belongs. 
Its explicit formula as well as statistical distribution was obtained and its rela- 
tion to the celebrated Riemann hypothesis was clarified [Zl El EJ ■ To prove the 
formula for fundamental cycle, one of the key steps is to compare a state with 
its 'reduced states' constructed by the '10-elimination'. In this article, we show 
that the initial value problem of the PBBS is solved by simple combinatorial 
arguments - essentially given in Ref. |Zj - with some remarkable features of the 
reduced states. 

First we quickly review the definition of the PBBS and its conserved quan- 
tities. Consider a one-dimensional array of boxes each with a capacity of one 
ball. A periodic boundary condition is imposed by assuming that the last box 
is adjacent to the first one. Let the number of boxes be N and that of balls be 
M. We assume M < N/2. An arrangement of M balls in N boxes is called a 
state of the PBBS. Denoting a vacant box by and a filled box by 1, a state 
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of the PBBS is represented by a 0, 1 sequence of length N. The time evolution 
rule from time step t to t + 1 can be described as follows: 

• For a given state, connect all 10 pairs in the sequence with arc lines. We 
call them 'l n arc lines'. 

• Neglecting the 10 pairs which are connected in the first step, connect all 
the remaining 10s with arc lines. We call them '2 n arc lines'. 

• Repeat the above procedure until all Is are connected to 0s with arc lines. 

• Exchange all the Is and 0s which are connected with arc lines. Then we 
obtain a new sequence which we call the state evolved by one time step. 

If we denote by Pj(t) the number of j n arc lines, we obtain a nonincreasing 
sequence of positive integers, Pj(t) (j = 1, 2, 3, . . . , to). Then, this sequence is 
conserved in time, that is, 

Pj{t) = Pitt + !) = Pj U = 1. 2, 3, ... , to). 

As the sequence (pi,P2, ■ ■ ■ ,Pm) is nonincreasing, we can associate a Young 
diagram to it by regarding pj as the number of squares in the jth column of the 
diagram. The lengths of the rows are also weakly decreasing positive integers. 
Let the distinct row lengths be L± > L2 > ■ ■ ■ > L s and let rij be the number 
of times that length Lj appears. The set {£j>^j}? = i is another expression for 
the conserved quantities of the PBBS. 

For example, for a state of the PBBS with N = 32, M = 14 

(ft) 00111011100100011110001101000000, 

the arc lines are drawn as in Fig. ^ an d its conserved quantities are expressed 
by the Young diagram given in Fig. [21 
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Figure 1: Arc lines and conserved quantities for the state (U). 



To solve the initial value problem of the PBBS, we utilize the fact that a 
state of the PBBS is determined by its conserved quantities and the positions 
of 10 pairs to which the T0-elimination' is applied. In the present context, a 
10- elimination is to convert a state to a state with smaller number of entries 
by eliminating all 10 pairs in the sequence connected with l n arc lines. Let 
us consider a state S with conserved quantities (pi,P2, ■ • • ,Pm) or equivalently 
{Lj,rij}j =1 . Note that s < to and L\ = m. We define its fc-reduced state 
[k = 1, 2, . . . , to) as the one obtained from the state by eliminating all the 10 
pairs connected with j n arc lines for all j : 1 < j < k. The length of the 
A:-reduced state is N — 2Y^j=iPj- We sometimes call the original state S as 
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Figure 2: Young diagram corresponding to the conserved quantities of (jj). 



0-rcduced state. Apparently 10-elimination is not a reversible operation. If, 
however, we remember the places where 10 pairs have been eliminated, we can 
recover the original state, at least up to shift, by inserting 10 pairs there. Since 
there has necessarily been a 10 pair between consecutive 1 and in the reduced 
state, we have only to remember the places for the other 10 pairs. (Such places 
are called the positions of 0-solitons in Ref. , and we use the same terminology 
here.) To specify the positions for insertion of 10 pairs, we number a place 
between j and j + 1th consecutive entries in a state by integer j. Note that, due 
to the periodic boundary condition, we set = N — 2 Y^j=i Pj f° r the fc-reduced 

j 

state. Hereafter we explicitly write the jth place by ' | ' in the 0, 1 sequence. 
For example, the state (jj) is expressed as 

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 

|0|0|1|1|1|0|1|1|1|0|0|1|0|0|0|1|1|1|1|0|0|0|1|1|0|1|0|0|0|0|0|0| . 

If we denote by E a 10-elimination, for the above state (jj), the fc-reduced 
state is expressed as -E fc (jj) and we have 

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 

£(jj) = |0|0|1|1|1|1|0|0|0|1|1|1|0|0|1|0|0|0|0|0|, 

1 2 3 4 5 6 7 8 9 10 11 12 13 14 

E 2 % = |0|0|1|1|1|0|0|1|1|0|0|0|0|0| , 

0123456789 10 

£ 3 (B) = |0|0|1|1|0|1|0|0|0|0| , 

1 2 3 4 5 6 

£ 4 (tf) = |0|0|1|0|0|0|, 

12 3 4 

£ 5 (») = 1 1 1 1 1 . 

In i?(jj) the positions of 0-soliton are 4, 7 and 15, in E 2 {$) the position of 
0-soliton is 10, in -E 3 (jj) there are no 0-solitons, and so on. Note that there 
is indeterminacy in constructing the original state from the reduced sequence; 
for example, both 0110001 and 1100010 turn to the same reduced state with 
the same positions of 0-solitons by 10-elimination. Hence we cannot necessarily 



determine the exact position of the 10 pair in the original state from the position 
of corresponding 0-solitons in the reduced state. 

In the fc-reduced states, E k (S), the number of positions of 0-solitons is pk — 

Pk+i (Pm+i '■— 0). Hence there appear V^^j = Pi 0-solitons in total in the 

i=i 

reduced states. We also find that 0-solitons appear only at the L^-reduced states 
E Lj (S) (j = 1,2, ... ,s) and that the number of 0-solitons is nj in E Lj (S). There 
is no 0-soliton in the other reduced states. So we denote by (k = 1,2, ... ,rij) 
the position of the fcth 0-soliton in the ij-reduced state. Since E Ll (S) only 
consists of N — 2M 0s, we can reconstruct the original state up to some shift 
from {Xfr}j' =1 2Li- H) however, we know the original position in the 0-reduced 
state of one of the 0-solitons in the Li-reduced state, we can recover the original 
state, because we have only to shift the state obtained by successive insertion 
of 10 pairs so that it coincides the original position. 

For this purpose, it is more convenient to introduce a set of variables {ol^ 
and X s+ i. First we define the reference positions Xj (j = 1, 2, . . . , s + 1) which 
is defined from x± recursively as follows. 

• We denote by Xli-i the position inbetween 1 and in the (L\ — 1)- 
reduced state of the 10 pair which turns into the 0-soliton at Xl 1 in the 
ii-reduced state. 

• Similarly we denote by Xi x _2 the position inbetween 1 and in the (L\ — 
2)-reduced state of the 10 pair which turns into the position Xl x -i in the 
(Li — l)-rcduced state. If there are more than one 10 pairs inserted at the 
position, Xl 1 _2 is the left most position among them. 

• Repeat the above procedure and obtain Xk (fc = 0, 1, ... , L\). 

• X 2 := X Lj (j = 1, 2, . . . , s + 1, where L s+1 := 0). 

In the above example ())), Xl 1 = X 5 = 2, X 4 = 3, X 3 =4, X 2 = 5, Xi = 6 
and X Q = 9. Hence X x = 2, X 2 = 3, X 3 = 5, X 4 = 6 and X 5 = 9. In the 
terminology of Ref. |7], Xj is the position of one of the 'largest soliton' in the 
.Lj-reduce state which turns to a 0-soliton in the ii-reduced state. Hence X s+ i 
is 'the original position in the 0-reduced state' of the 0-soliton at the position 
r (i) 

Then we define a[ 3) (1 < a k j) < Nj) by 

4 j) =X j -x[ j) mo&Nj (i = l,2,...,s, k = l,2,..., nj ), (1) 

where 

j 

Nj :=N- 2M + 1n t (L t - Lj) (j = 1, 2, . . . , s) 

i=l 

is the number of entries in the L^-reduced state. Note that = N\ and 

is the distance of the fcth 0-soliton from the position of the 'largest soliton' in 

the Lj-reduce state, which turns to a 0-soliton in the Li-reduced state. Since 
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the state S can be determined up to shift by inserting 10 pairs at the posi- 
tions of 0-solitons and that between consecutive 1 and 0, and X s+ \ determines 
the amount of the shift, the state S is uniquely determined by the variables 
{oik ^Li+Si an d X s+ \. Formally we may write that 

} TkU e snj ( z *i ) := Ziv, x Zat. x • • • x Zat. (ST* (j = 2, 3, . . . , s), 

where Z ^ j is the cyclic group of order Nj , S rlj is the symmetric group of order 
rij, and, since there are iVi + 1 distinct positions 1 for (n\ — 1) 0-solitons in the 
Li-reduced state, 

and X s+ i s In- If we define 

V Y ■= S ni - l {1 Nl+ x) x 5 n2 (Zjv 2 ) x ■■• x S n *{Z Ns ) x Z N , 

then an element of Vy naturally corresponds to a state of the PBBS. Since 
there are n\ choices of the reference position, there are exactly n\ elements in 
Vy which correspond to a state of the PBBS. We regard that one element in Vy 
is equivalent with another if and only if they corresponds to the same state of 
the PBBS. Denoting by Vy the quotient set of Vy according to this equivalence 
relation, we obtain the following theorem; 

Theorem 1 

Denote by ily a set of the states of the PBBS of M balls and N boxes with 
conserved quantities {Lj 7 Uj}j =1 characterized by the Young diagram Y . Then 
there is a one to one correspondence between an element of fly &nd that ofVy, 
hence Qy — Vy. The explicit bijection is given by the 10-eliminations and its 
inverse operations with shift using the variables X s+ i and }j' =1 ^L 1+Sl . ■ 

As for the explicit construction of a state of the PBBS from an element in Vy, 
see the example given below. 

Now we consider the initial value problem of the PBBS. Let S(t) be the state 
evolved from the state S by t time steps. From Theorem^ we have that the dy- 
namics of the PBBS can be described by an element in Vy. Hence, to determine 
a state S(t), it is enough to obtain the variables {^{ct^ {t)Yj=\ ^i+Sx ■> ^s+l(X)) 

Vy from the initial values (^{a^\o)} S j' =1 ^ 1+Si X s+ i(0)j € Vy. However, the 
time dependence of these variables has already been given in Ref. [7]: 

Proposition 1 (0 Theorem 3.1, Lemma 4.2 and Lemma 4.3) 

For i = 1, 2, . . . , s + 1, let ^(t) be 71 (t) := 0, 72 (t) := (Li - L 2 )t and 

i-l rij 

7i(t) := (£1 - U)t + 2J2(Lj - Li) E 

J=2 fc=l 

^^When another 0-soliton is located at position Xi , it is either at the left of the reference 
0-soliton or to the right of it. 
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for i = 3, 4, . . . , s + 1, where 

7,-w+4 j) (o)-i 



N 3 



(j = 2,3,...,s, fc = l,2, 



Then it holds that 

4 J) W=4 J) (0)+7,W modiV,-. 

and 

X s+1 (t) = X s+1 (0) + 7s+ i (f ) mod N. 

Remark 1 

Note that both Ji(t)s and a^(t)s are determined recursively. As we shall see, 
practically we have only to use the relation 

4 i) (o)+7iW = iv J 4 j) w+4 i) (t). 

In conclusion, we have solved the initial value problem of the PBBS which 
may be stated as 

Theorem 2 

The initial value problem of the PBBS is solved in the space ofVy ■ Its dynamics 
is explicitly given in Proposition^ 

In the rest of this article, we explain how to obtain the time evolution of a 
state of the PBBS by means of an example. Suppose that we have the state (|) 
which we have used previously 

(0) 00111011100100011110001101000000 

at time step t = 0. From Fig.^and £ fc ()}), we find s = 4, Ni = 4, N 2 = 6, iV 3 = 
14, AT 4 = 20. Then referring to the reduced sequences E k (§), we obtain 

x[ 1] (0) = 2; 4 2) (0) = 3; xf\o) = 10; x^ } (0) = 4, 4 4) (0) = 7, 4 4) (°) = 15 - 
The reference positions Xj(0) (j — 1,2, ... , 5) are found to be 

Xi(0) = 2, X 2 (0) = 3, X 3 (0) = 5, X 4 (0) = 6, X 5 (0) = 9. 
Hence a^ } (0)s are given as 

«P(0) = 4; a^ 2) (0) = 6; a< 3) (0) = 9; a^{0) = 2, (0) = 19, 4 4) (0) = H. 
Now let us consider the state at t = 10000. According to Proposition^] oc£\t)s 
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and 7,(t)s are calculated recursively as follows; 

(1) N 2 = 6, 72 (t) = (Li -L 2 )t= 10000 

=► af } (0) + 72 0) = 6 + 10000= 1667-6 + 4 
=* af } (i)=4; 

(2) iV 3 = 14, 73 (i) = (Li - L 3 ) • t + 2{L 2 - L 3 ) • 1667 = 36668 

=► af } (0)+7 3 (i) = 9 + 36668 = 2619-14+ 11 
af ) (i) = ll; 

(3) N A = 20, 

74 (i) = (ii - L 4 ) • t + 2{L 2 - Li) ■ 1667 + 2{L 3 - L A ) ■ 2619 = 55240 
of (0) + 74 (i) = 2 + 55240 = 2762 -20 + 2 
a 2 4) (0) + 74 (t) = 19 + 55240 = 2762 -20 + 19 
4 4) (°) + H(t) = 11 + 55240 = 2762 -20 + 11 

ctf\t) = 2, 
a 2 4) (*) = 19, 



Q 



(*) = ii; 



(4) N = 32, 

7s(t) = (Li - L B ) • t + 2(L 2 - L B ) • 1667 + 2(L 3 - L 5 ) ■ 2619 
+2(L 4 - L 5 ) ■ (2762 + 2762 + 2762) = 90384 
=> X 5 (t) = X 5 (0)+~/ 5 (t) mod N = 9 + 90384 mod 32 
= 25. 

From these data, the state at t = 10000 is constructed up to shift by inserting 
10 pairs as; 

1 2 3° 4 

1 1 1 | 1 , 

0° 1 2 3 4* 5 6 

— |0|0|0|1|0|0|, 

0123456 7* 89 10 

— > |1|0|0|0|0|1|1 1 1 1 | , 

1 2 3 4 5 6 7 8 9 10* 11 12 13° 14 

|1|1|0|0|0|0|0|1|1|1 | 0|0|0 1 0| , 

01 2° 3456789 10 11° 12 13* 14° 15 16 17 18 19 20 

— > 1 1 1 1 | 1|0|0|0|0[0|0|1|1 | 1 1 1 | | 0|0|0|1|0|0| , 

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20* 21 22 23 24 25 26 27 28 29 30 31 32 

— > |1|1|1|0|1|1|0|0|0|0|0|0|0|1|1|1|0|1|1|1|0|0|1|0|0|0|0|1|1|0|0|0| . 

3* 

Here ' | ' denotes the position of the largest soliton (referring position) and 

' | ' denotes that of a 0-soliton. The position x\ (t) can be chosen arbitrary, 

and we took x^p(t) = 3 in this example. Finally we translate the above state 
so that the position of the largest soliton coincides with X s+1 (t) as 

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 34 25* 26 27 28 29 30 31 32 

|1|1|0|0|0|1|1|1|0|1|1|0|0|0|0|0|0|0|1|1|1|0|1|1|1|0|0|1|0|0|0|0| . 
This is the state at t = 10000. 
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In this article, we have solved the initial value problem in an elementary 
way. We also remark that our method can be equally applied to extended 
PBBSs with carrier capacity £ as those treated in Ref. \T2\. I n these PBBSs, 
we have only to replace Lj — > miii[Lj,£] in Proposition ^ and apply the above 
procedures. Clarifying the relation between our methods and previous work 
based on algebraic curves and representation theories, is one of the important 
problems we want to address in the future. 
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